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At low temperature, a thermodynamic system undergoes a phase transition when a physical parameter 
passes through a singularity point of the free energy. This corresponds to the formation of a new order. At 
high temperature, thermal fluctuations destroy the order. Correspondingly, the free energy is a smooth 
function of the physical parameter and singularities only occur at complex values of the parameter. Since a 
complex valued parameter is unphysical, no phase transitions are expected when the physical parameter is 
varied. Here we show that the quantum evolution of a system, initially in thermal equilibrium and driven by 
a designed interaction, is equivalent to the partition function of a complex parameter. Therefore, we can 
access the complex singularity points of thermodynamic functions and observe phase transitions even at 
high temperature. We further show that such phase transitions in the complex plane are related to 
topological properties of the renormalization group flows of the complex parameters. This result makes it 
possible to study thermodynamics in the complex plane of physical parameters. 



T 



he physical properties of a many-body system in the thermodynamic equilibrium are fully determined by the 
partition function Z(j5,/li,/l2,. . which is a function of the coupling parameters {1^} of the system and the 
temperature T (or the inverse temperature = 1/T). The partition function is the summation of 



the Boltzmann factor e 



= Trk 



where the 



over all energy eigen states, i.e., Z = 

Hamiltonian H = k^k^k is characterized by a set of coupling parameters {kj} (e.g., in spin models the 
magnetic field h = Xi, the nearest neighbor coupling / = /I2, and the next nearest neighbor coupling /' = >^3, 
etc.) . The normalized Boltzmann factor Z~ ^exp( — PH) is the probability that the system is in a state with energy H. 

At a sufficiently low temperature (P > \^k\~'^)y the many-body system may have different orders (phases) in 
different parameter ranges and therefore phase transitions may occur when these parameters are varied (see 
Fig. la). Phase transitions correspond to non- analytic or singularity points of the free energy F (which is related to 
the partition function by F = —p~^\nZ). For a finite-size system, the partition function is always positive and 
hence the free energy always analytic for real physical parameters. The partition function can be zero for complex 
temperature or complex coupling parameters (see Fig. lb & Ic). Such is the case for Lee- Yang zeros^'^ in the 
complex plane of magnetic field for a spin lattice and Fisher zeros^ in the complex plane of temperature. At a 
sufficiently low temperature and in the thermodynamic limit (where the number of particles in the system 
approaches infinity), the zeros of the partition function can approach to real axes of the coupling parameters, 
and, hence, phase transitions can occur ^'^ (see Fig. Ic). At high temperature, the free energy would be non-analytic 
only for complex values of temperature and coupling parameters'"^ (see Fig. lb). Therefore, no phase transitions 
would occur (see Fig. la). Physically, this is because, at high temperatures, thermal fluctuations destroy all 
possible orders and prevents the phase transitions. Phase transitions at high temperatures would be possible if 
one could access complex parameters, which, however, are generally regarded as unphysical. 

Here we present a systematic approach to the complexification of an arbitrary physical parameter, using a 
quantum evolution that involves a complex phase factor. There are two fundamental laws of physics that involve 
probabilistic distributions. One is the Boltzmann distribution Z~^exp{ — fiH), the probability of a system being in a 
state with energy H at inverse temperature p. The other is the Schrodinger equation that gives exp( — /YH), the 
probability amplitude of a system being in a state with energy H at time t. The unit imaginary number / = \/— 1 is 
associated with the time in quantum evolution. Thus a time-dependent measurement of the system 



L(0=Z-^Tr[exp(-j5H-/YH/ 



(1) 



has the form of a partition function with an imaginary coupling parameter itip. By choosing Hj to be the k-\h 
component ofH = ^ k^kHh i-^o Hj = Hj^, the k-th coupling parameter is analytically continued to the complex 
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Figure 1 | Phase transitions and zeros of a partition function in the 
complex plane of a physical parameter, (a). The schematic phase diagram 
of a general many-body system. At low temperature, the system has two 
ordered phases, I and II, separated by a critical parameter 2^- At high 
temperature, the system is in the disordered phase. When we tune the 
control parameter at low temperatures (as indicated by blue dashed arrow) 
we cross two phase boundaries (indicated by solid-red lines) and therefore 
experience two phase transitions. Conversely at high temperature 
sweeping the parameter in the same range (blue dotted arrow) would not 
cause a phase transition, (b). The schematic distribution of zeros of the 
partition function in the complex plane of the parameter for temperature T 
above the critical point T^. The two blue points mark the edge singularities, 
(c). The same as (b) but for the temperature below the critical point ( T < 
Tc). The singularity edges approach to the real axis in the thermodynamic 
limit. 



plane via ^k~^ it/p. One can also realize complexification of the 
inverse temperature (P ^ + it) by choosing Hj = H. Such a time- 
domain measurement, as will be discussed later in this paper, is 
experimentally realizable, though non-trivial in general. The time- 
domain measurement in equation (1) is fully determined by the 
equilibrium-state partition function of the system. The measurement 
result will be zero when a zero of the partition function in the com- 
plex plane is encountered^. Lee-Yang zeoros have been recently 
observed via such measurement^. Furthermore, critical times may 
exist corresponding to the singularity points in the complex plane, 
as indicated by recent studies of central spin decoherence caused 
by an Ising bath^ and evolution of a quantum spin modeP. The 
measurement in equation (1) provides a systematic approach 
to accessing different types of zeros in the complex plane for the 
partition function, such as Lee-Yang zeros^'^ by choosing 
Hi = ^ . for a system of spins {S^j, Fisher zeros^ by choosing Hj 
= H^, and other types of zeros (yet unnamed) by choosing 
Hi = ^ ^J^J+v ^ ^ ^j'^j + 1' ^t^- More importantly, by com- 
plexification of a physical parameter, we have a way to access the 
singularity points in the complex plane and to reveal phase transi- 
tions at high temperature (j5"^ > {\^k\}) via time-domain measure- 
ment. The high-temperature phase transitions are also inferred from 
a recent finding that quantum criticality can emerge at high tem- 
peratures by long-time quantum evolution^. Note that the phase 
transitions studied in this paper are not the time- domain or dynam- 
ical phase transitions studied in Refs. 6 & 8. Instead, they occur 
between different parameter regimes, signified by the qualitative 
changes of time-domain measurement. 



High-temperature magnetic phase transitions. Spin systems can 
have ferromagnetic (FM) or antiferromagnetic (AFM) orders at low 
temperatures, corresponding to positive or negative coupling (/) 
between the spins, respectively. Thus an FM-AFM transition 
would occur if the coupling / varies from positive to negative. At 
high temperatures, thermal fluctuations destroy the magnetic order 
and hence no phase transition is expected with changing /. Here we 
study the Ising spin model to demonstrate the FM-AFM transition in 
the complex plane of the parameter /. The Hamiltonian for the 
general Ising model is 

where is the coupling between spins ct/ and Gp h is the magnetic 
field, and the spins Gj take values ± 1. At low temperatures, when the 
coupling changes from positive to negative, the Ising model presents 
a phase transition from the FM order to the AFM order at zero field. 
Correspondingly, the Lee-Yang zeros in the complex plane of scaled 
magnetic field z = exp(2ph) exhibit different distributions. Note that 
the distinct features of Lee-Yang zeros distribution in the complex 
plane persist even at high temperatures (T > |/y|). 

To be specific, we study the one- dimensional (ID) Ising model 
with nearest-neighbor coupling /, which can be exactly solved 
through the transfer matrix method^°"^^ (see Supplementary 
Information). There is no finite temperature phase transition in 
the ID Ising model. The Lee- Yang zeros of the ID Ising model of 
N spins have been exactly calculated^. We plot the distribution of 
Lee- Yang zeros in Fig. 2a. For AFM coupling (/ < 0), all the zeros lie 
on the negative real axis (indicated by the red surface) (Fig. 2a). 
While for the FM coupling (/ > 0), the zeros are distributed on an 
arc of the unit circle (indicated by the blue surface) (Fig. 2a). At the 
transition point (/ = 0), all the Lee- Yang zeros are degenerate at z„ = 
— 1 (indicated by the green solid ball). 

To observe the Lee-Yang zeros and the critical behaviors in the 
complex plane^'^, we study the time- domain measurement in equa- 
tion (1) by choosing Hj = ^j^j ^^d analytically continuing the 
magnetic field to the complex plane h + it/ p. In this case, the 
time-domain measurement corresponds to decoherence of a 
quantum probe spin coupled to the Ising modeP'^ (see physical rea- 
lizations). Figs. 2b-e plot the time-domain measurement calculated 
for different coupling parameters / at fixed inverse temperature P = I 
(which is a high temperature case for the ID Ising model, since the 
critical temperature of this model is zero). For AFM coupling / = — 1, 
the time-domain measurement has no zeros (Fig. 2b) and it is a 
smooth function of time. On the contrary, for the FM coupling / = 
1, the time-domain measurement shows a number of zeros (Fig. 2c), 
which have a one-to-one correspondence to the Lee-Yang zeros^. 
Approaching the thermodynamic limit, the Yang-Lee edge singular- 
ities (the starting and ending Lee- Yang zeros along the arcs)^'^ lead to 
critical times in the time-domain observation^. To demonstrate this, 
we perform a finite size scaling analysis on the time-domain mea- 
surement and show the scaled results \L{t)\^'^ in Figs. 2d & 2e. The 
profiles of the time-domain measurement in the FM and AFM 
regions are qualitatively different. For AFM coupling, the scaled 
measurement is a smooth function of time (Fig. 2d). While for FM 
coupling, the scaled measurement presents sudden changes at critical 
times corresponding to the Yang-Lee edge singularities (Fig. 2e). The 
profiles of the time-domain measurement in the FM and AFM 
regions cannot be smoothly transformed into each other, which is 
the signature of the onset of a high-temperature phase transition. 

We further study the high-temperature AFM-FM phase transition 
in the time-domain measurement for a two-dimensional (2D) Ising 
model. Specifically, we consider a 2D Ising model in a square lattice 
with nearest neighbor coupling /. This model under zero field is 
exactly solvable^ ^'^^ and has a finite-temperature phase transition at 
Pc ~ 0.44/ Fig. 3 shows the time-domain measurement in the 2D 
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Figure 2 | High-temperature AFM-FM phase transitions revealed by 
time-domain observation, (a). Schematic 3D plot of the Lee-Yang zeros in 
the ID Ising model. For / > 0, the Lee-Yang zeros are distributed on an arc 
of the unit circle with two edge singularity points determined by the 
coupling strength (blue surface). As the coupling strength decreases the 
arcs shrink and finally the singularity edges merge to a point at z = — 1 for / 
= 0 (indicated by the green solid ball). For / < 0, the Lee-Yang zeros lie on 
the negative real axis (red surface), (b). The time-domain measurement as 
a function of time for ID Ising model with N = 20 spins at P = 1 and / = 
— 1 (AFM). (c). The same as (b) but for / = 1 (FM). (d). Finite size scaling 
of the time-domain measurement in the ID Ising model at P = 1 and / = 
- 1 (AFM), the red Une is for AT = 100 spins, green line for N = 500 and 
blue line for N = 1000; (e). The same as (d) but for / = 1 (FM). 

Ising model calculated for different coupling parameters / at a fixed 
high temperature, = 03 < ^c- For the AFM coupling / = — 1, the 
partition function has no Lee- Yang zeros on the unit circle {\z\ = 
|exp(2j5/z)| = 1)^^ and therefore the time-domain measurement has 
no zeros (Fig. 3a). While for the FM coupling / = 1, the time-domain 
measurement presents a number of zeros (equal to the number of 
spins) (Fig. 3b), corresponding to the Lee- Yang zeros along the unit 
circle. We do finite-size scaling analysis in Fig. 3c and Fig. 3d for 
AFM coupling and FM coupling, respectively. It is clear that the 
time-domain measurement is a smooth function of time for the 
AFM coupling, while it presents sudden changes at critical times 
corresponding to the Yang- Lee edge singularities for the FM coup- 
ling (Fig. 3d). Thus a phase transition with varying the coupling 
constant / occurs at a temperature higher than the critical temper- 
ature (j5 < 0.44/|/|). 

Renormalization group theoretic analysis. The renormalization 
group (RG) method, a powerful tool for studying conventional 
phase transitions, can be applied to the phase transitions in the 
complex plane of physical parameters. Since the phase of any 
complex number is defined modulo 271, the RG flows of complex 
parameters can present novel topological structures. 

As an example we first consider the ID Ising model and define the 
dimensionless parameters iCo = j^/and/zo — Ph.The renormalization 
of the model can be exactly formulated by blocking two neighboring 




Figure 3 | High-temperature AFM-FM phase transitions of the 2D Ising 
spin model in a square lattice revealed by time-domain observation, a. 

The time-domain measurement as a function of time for the 2D Ising 
model in a AT = 4 X 5 square lattice, ^ = 0.3 and / = - 1 (AFM). b. The 
same as (a) but for / = 1 (FM). c. Finite size scaling of the time-domain 
measurement with fi = 0.3 and / = — 1 (AFM), the green line is for AT = 4 X 
100 spins, the blue line is for AT = 6 X 100 and red line for AT = 8 X 100. d. 
The same as (c) but for / = 1 (FM). 

spins into one (Fig. 4a) By continuation of the dimensionless 
external field to a purely imaginary value ho = /to/2 (the physical 
external field is zero), the exact RG flow equations become^^ 



Ti(To,Ko) = To- 



. / cosh(2K'o + ho / 2) 
^ ^ Vcosh(2K'o-/To/2) 



I f cosh(4iCo) + cos To 

Ki(to,Ko)= -In — 

4 \ 1 + cos To 



(3) 



where the coupling Ki remains real after renormalization. Since t is 
defined modulo In, the parameter space can be identified with the 
surface of an infinitely long cylinder with unit radius. The original 
system corresponds to the parameter curve K= Ko and — 7i < Tq < 7i, 
and therefore its winding number ( W#), defined to be the number of 
times the parameter curve wraps around the cylinder, is 1 (Fig. 4b). 

The parameter curve is renormalized according to the RG flow 
equations. Since Ki(to,Ko) = Ki(to,— iCo) in Eq. (3), the distinct 
behaviors of the FM and AFM Ising chains are encoded entirely in 
different RG flows of t. This is illustrated in Figs. 4c-h. Figs. 4c, 4e & 
4g show that after successive renormalization, the winding number 
in the FM case (Ko = 1/8) becomes 2, 4, 8,.... On the contrary, the 
winding number in the AFM case (Ko = — 1/8) is zero after renor- 
malization (Figs. 4d, 4f Sc 4h). The RG flow equations become trivial 
at the phase transition point Ko=0, which corresponds to the infinite 
temperature limit, and the winding number remains unchanged 
(W!f^ = 1) after the renormalization. In summary, the winding num- 
bers for different couplings after k steps of renormalization are 



0 for Ko<0 

1 for Ko = 0. 
2^ for Ko>0 



(4) 



The different topologies of the RG flows demonstrate unambigu- 
ously the high-temperature phase transition with varying the coup- 
ling parameter. 

We further consider the RG flows of the 2D Ising model in a square 
lattice. By continuation of the external field to a purely imaginary 
value of h = /t/2, the approximate RG flow equations read^^ (See 
Supplementary Information for derivation) 
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Figure 4 | Renormalized parameters of the ID Ising model with a purely imaginary field, (a), Real-space renormalization scheme. In each 
renormalization step the spins on every other site are traced out, which effectively combines pairs of neighboring spins into block spins with renormalized 
coupling and external field, (b), The un-renormalized parameters. The spin coupling iCand evolution time t (i.e., imaginary part of the external field) 
form a space identical to the surface of an infinitely long cylinder. When the imaginary field t is varied at a fixed value of iC = Xq, the curve winds about the 
cylinder once and so the winding number W# = 1. (c-h). Under RG flow, the original curve depicted in (a) would transform differently in the different 
parameter regimes of Xq. (c), (e) and(g) in turn show the parameters for the FM case with iCo = 1/8 after one, two, and three steps of renormalization. (d), 
(f) and (h) show the corresponding results for the AFM case with Xq = — 1/8. The vertical dashed lines indicate t = 0 mod (2 71), which are identified as the 
same line when represented on an infinitely long cylinder. The winding numbers can be directly inferred from the number of times the renormalized curve 
crosses this line. 



r , i, /cosh (4iCo + /To/2) cosh' (2Ko + /To /2)\ 

Ti(To,iCo)~To m ^ 

^ 2 Vcosh(4Ko-/To/2)cosh'(2Ko-/To/2)y 

I ^ ^ 16 V 1+ cos (To) J 

Fig. 5 presents the RG flows of the parameters in the FM and AFM 
cases. Figs. 5a, c & e present the renormalized parameters under one, 
two and three applications of the renormalization transformation for 
the FM case (Kq =1/8). Figs. 5b, d & f present the results for the AFM 
case (Kq = — 1/8). The winding numbers of the different cases after k 
steps of renormalization are 

r [(5/2/] for Ko>0 
W^^l 1 for Ko = 0, (6) 

[ [(5/2/-V2] for Ko<0 

where [x] is the integer part of x. The winding number of the para- 
meters under RG reflects the different topologies intrinsic to the RG 
flow equations in the different parameter regimes. 

Transverse-field Ising model. The models considered above are all 
classical models in which different components of the Hamiltonian 
commute. A natural question arises about whether the high- 
temperature phase transitions with complex parameters would 



exist also for quantum models. To address this question, we study 
the ID transverse-field Ising model. The model contains N spin- 1/2 
with nearest neighbor interaction (li) along the x-axis and under a 
transverse field {X2) along the z-axis, described by the Hamiltonian 

N N 

H = hH, + = hY, ^"^iV 1 + ^2 ^ ^tJ, (7) 

where a^J^^^ is the Pauli matrix of the j-th spin along the x/y/z-axis. 
This model is exactly solvable^^. It has a quantum phase transition 
between a magnetic ordered phase for \ki\ > I/I2I and a disordered 
phase for < |>^2| at zero temperature, but has no finite- 
temperature phase transition for any parameters on the real axis. 
By defining the dimensionless magnetic field h = ^2^1iy the Lee- 
Yang zeros are determined by Re(hy + lm{hy =1 + [{n + 1/2)71/ 
PY, |Re(/z)| ^ 1. Therefore, the zeros are located on circles and have 
cutoff at singularity edges Re(/z) = ±l,lm(h) = ±\(n + l/2)n/p\ (see 
Fig. 6a). When the temperature approaches zero (|j5| ^ 00), the Lee- 
Yang zeros are on the unit circle. When the temperature is high 

< 1), the radii of the circles are ^1 + + l/2)7r/j5|^ » 1. 
Therefore, the zeros are distributed, approximately, along 
horizontal lines with interval \n/p\. The fact that the Lee- Yang 
zeros exist only in the parameter range of\h\ < 1 indicates that the 
time-domain measurement of the system would present phase 
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Figure 5 | Topological difference between the renormalization group 
flows of the FM and AFM 2D Ising models, (a), (c) & (e) present in turn 
the parameters (spin coupling and imaginary external field) after one, two, 
and three applications of the renormalization transformation in the FM 
case {Kq = 1/8). (b), (d) & (f) are the same as (a), (c) & (e) in turn but for 
the AFM case (Xq = — 1/8). The renormalized parameters for the FM and 
AFM cases have different winding numbers (determined by how many 
periods (2%) of time the parameter curves are spanned over). 

transitions between the two parameter regions, |>^2| > I'^il and I/I2I < 
|>^i I . Specifically, the time-dependent measurement can be devised as 
L{t) = Z"^Tr[exp( — /tH2)], which is the partition function with 
a complex external field (I2 ^2 ~^ it/P). The contour plot of the 
time-domain measurement as a function of external field and time 
are presented in Fig. 6b. To demonstrate the phase transitions more 
clearly, we plot in Fig. 6c the time- domain measurement as a 
function of external field for different times. Since the zeros are 
bounded in the range \h\ ^ 1 with Yang-Lee edge singularities at 
Re{h) = ±1, lm{h) = ±\{n + 1/2)71/ the time-domain 
measurement has a sharp change when we tune the parameter 
from \h\ < 1 to |/z| > 1 at times lm{h) = ±\{n + \l2)%ip\ (see 
Fig. 6c). 

Physical realization. The time-domain measurement in equation (1) 
resembles the Loschmidt echo^^'^^, or equivalently, decoherence of a 
central spin coupled to the system. Thus we may implement the time- 
domain measurement by coupling the system to a central spin 
through the probe-system interaction 1 1> (T | ® ^^t(0 + I i) (i I ® 
i'^x = I TXT I ~ ai^d measuring the central spin coherence. 

Essentially, the coherence of the probe spin is a complex phase 
factor associated with a real Boltzmann probability for each state of 
the system. Therefore the probe spin coherence measurement 
amounts to continuation of a physical parameter to the complex 
plane. If we initialize the central spin in a superposition state | T) + 
II) and the system in a thermal equilibrium state described by the 
canonical density matrix, p = Z"^Tr[exp(— jSH)], the probe spin 
coherence is 



-20-^ , ^ , ^ , ^ . , 

-2-1012 

Re(/z) 




Figure 6 | Finite-temperature phase transitions of a transverse field Ising 
spin chain revealed by time-domain observation. The one- dimensional 
model contains 100 spins. The temperature is such that P = 0.5. (a). Lee- 
Yang zeros in the complex plane of the external field, (b). Contour plot of 
the time-domain measurement, IglLl, as a function of the external field 
Re(h) and the time lm(h). (c). IglLl plotted as functions of the external 
field Re(/z) for different times lm{h) = 3.10 (red-dashed line), 3.14 (black- 
solid line), 3.26 (blue-dashed line) and 3.35 (green-dashed line). The phase 
transitions at Re(/z) = ±1 are observed. 



(S+ ) = (S,) + /(S^) =Tr[e-^«e'-(«+«T)'e-'(H+Hi)t] /Tr[e-/'«] . (8) 

If [HiH] = 0, the time-domain measurement L{t) = 
Z"^Tr[exp(— j5H)exp( — /H/)], which is equal to the probe spin 
coherence with a probe-bath coupling Hsb = Sz(S)Hi, i.e., H| = 
-H| = Hi/2. Or if [HjH] 7^ 0 but [[HiH],H] = [[HjHl Hj] = 0, 
the time-domain measurement can also be factored as L(t) = 
Z"^Tr[exp(— j5H)exp( — /Hjt)] (See Supplementary Information for 
details) and can be implemented by the probe spin coherence with 



SCIENTIFIC REPORT: | 4 : 5202 | DOI: 1 0.1 038/srep05202 



5 



a modified probe-bath coupling = Hj/2 — H &:H^ = —Hj/2 — H. 
If [H,Hi] = isHj, the time-domain measurement can be written as L(t) 
= Z"^Tr[exp(— sm(s)Hj/s)exp(—pH)] (See Supplementary 
Information for details), which can be implemented by probe spin 
decoherence with a probe-bath coupling = sin(s)Hj/2s — H 
= —sm{s)Hj/2s — H. In the general case, the time-domain mea- 
surement in equation (1) can be written as L{t) =Z~^TY[e~^^ 

e-'^'^'^l with e-'^^ = Texp(^-^|'H(w)^lw^Texp^-0'H(^ 

exp(- itH'i)= exp(j5H')exp(- PH-itHi) and H{u) = exp{iutHi/2)H 
exp{ — iutHi/2), where T and f are the time-ordering and anti- 
ordering operators, respectively. If one initializes the bath in a 
canonical state p = exp( — j5H')/Tr[exp( — the time-domain 
measurement can be implemented by probe spin decoherence with 
probe-bath coupling H| =H'i/2-H' kH^ = -H'/2-H', up to a 
normalization factor (see Supplementary Information for details). 
The physical realization of the modified Hamiltonians H' and H'l 
is non-trivial. 

Note that the time-domain measurement in equation (1) is similar 
to the measurement of the characteristic function of the work distri- 
bution in a quantum quench^^'^\ which plays a central role in the 
fluctuation relations in non- equilibrium thermodynamics^^"^^. The 
probe spin decoherence realization of the time-domain measure- 
ment can also be related to the quench dynamics where the evolution 
of the system can be controlled under different Hamiltonians for 
different periods of time^^. 

Summary. We have shown that the quantum evolution of a system 
originally in thermodynamic equilibrium is equivalent to the 
partition function of the system with a complex parameter. By 
choosing different forms of coupling we have a systematic way to 
realize the continuation of an arbitrary physical parameter to the 
complex plane. The time-domain measurement allows us to study 
a variety of zeros of the partition function. More importantly, we can 
access the singularity points of thermodynamic functions in the 
complex plane of physical parameters and therefore observe phase 
transitions at high temperatures. The physical realization of the time- 
domain measurement may be nontrivial but in principle it can be 
implemented by probe spin decoherence or quantum quench 
experiments. This discovery makes it possible to study 
thermodynamics in the complex plane of physical parameters. 

Methods 

The ID Ising spin model was exactly diagonalized by the transfer matrix method^°'^\ 
The evaluation of the partition function after the transformation becomes a trivial 
problem of diagonalization of a 2 X 2 matrix. The probe spin coherence (which has 
been formulated in terms of partition functions^) was similarly calculated. Similarly, 
by applying the transfer matrix method, the 2D Ising model without magnetic field 
was mapped to a ID Ising model with a transverse field^°'^\ For a 2D Ising model in a 
finite magnetic field, it was mapped by the transfer matrix method to a ID Ising model 
with both longitudinal field and transverse field^\ which was then numerically 
diagonalized. Therefore, the partition function and hence the probe spin coherence 
for the 2D Ising model in finite magnetic field were obtained. We derived the exact RG 
equations of the ID Ising model and approximate RG equations for the 2D Ising 
model in square lattice for real parameters^^'^^. By analytic continuation we obtained 
the RG equations in the complex plane of the physical parameters and analyzed the 
RG flows of the complex parameters. The ID Ising model with a transverse magnetic 
field was exactly solved^*^ and the partition function and time-domain measurement 
were then calculated. 

Full methods are included in the Supplementary Information. 
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